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We study the large deviations statistics of the intensive work done by changing globally a control 
parameter in a thermally isolated quantum many-body system. We show that, upon approaching 
a critical point, large deviations well below the mean work display universal features related to the 
critical Casimir effect in the corresponding classical system. Large deviations well above the mean 
are, instead, of quantum nature and not captured by the quantum-to-classical correspondence. For 
a bosonic system we show that in this latter regime a transition from exponential to power-law 
statistics, analogous to the equilibrium Bose-Einstein condensation, may occur depending on the 
parameters of the quench and on the spatial dimensionality. 



PACS numbers: 

Introduction - Recent experimental progresses in the 
physics of trapped ultracold atomic gases have stimu- 
lated a growing interest in the non-equilibrium behavior 
of thermally isolated quantum many-body systems [1 , . A 
number of aspects are presently being investigated exper- 
imentally, ranging from the propagation of correlations 
after quenches ^ to relaxation and pre-thermalization 
inferred from the statistical fluctuations of the interfer- 
ence contrast of split condensates [J . On the theoreti- 
cal side, a compelling issue under investigation is that of 
the role played by universality in the non-equilibrium dy- 
namics ^ , since predictions independent of microscopic 
details make the comparison with experiments a partic- 
ularly stringent test. Universal behavior can be investi- 
gated by studying either the time dependence of correla- 
tion functions [H [S] , in particular close to criticality, or 
their statistical fluctuations [3]. In this context, a number 
of studies have focused on macroscopic, thermodynamic 
variables such as work |7H12j and entropy jl3j , exploring 
the emergence of universality in their statistical fluctua- 
tions. 

Statistical fluctuations are known to provide in- 
sight into the physics of classical equilibrium and non- 
equilibrium systems [13]. The statistics of macroscopic 
extensive variables exhibits a first, obvious form of 
universality associated to typical, "small" fluctuations, 
which is however rather insensitive to the underlying 
properties of the system [T?]. Indeed, as the mean of 
a generic extensive quantity Wn (e.g., the magnetization 
in a spin system) grows proportionally to the number N 
of degrees of freedom, the one of the associated inten- 
sive variable wn = Wn/N (i.e., the magnetization per 
unit volume) approaches a flnite value w. The central 
limit theorem (when applicable) suggests that the typi- 
cal fluctuations of wn are suppressed ~ l/N^/"^ and have 
a Gaussian distribution around w. On the other hand, 
large fluctuations, though rare, are capable of probing 
the specific details of the physical system and they 



might provide valuable information on its universal be- 
havior. In order for a large fluctuation to occur, an ex- 
tensively large number of microscopic fluctuating vari- 
ables (i.e., the spin, in our example) has to deviate sig- 
niflcantly from their corresponding means [TS] and this 
happens with a probability which is exponentially small 
in the size N . Accordingly, for large N , one expects 
Wn to be distributed according to a probability density 
p{w) ~ eyLp[—NI[w)], where the so-called rate function 
I{w) is non-negative, vanishes for w = w |20| . and char- 
acterizes the statistics of both large deviations and Gaus- 
sian fluctuations. 

Here we show that the statistics of large deviations of 
the intensive work w done during a global quench of a 
thermally isolated quantum many-body system provides 
insight into its universal properties. For a global quench 
one heuristically expects p{w) to feature a prominent 
Gaussian peak centered at a finite mean w. By focus- 
ing on the tails of this distribution, we demonstrate that 
there is a clear distinction between large deviations well 
below (w <C w) and well above {w ^ w) the mean. The 
former are determined by the excess free energy f^x of the 
d ~\- 1 dimensional classical correspondent in a film jB] |5] 
and acquire universal features close to a possible critical 
point. The latter, instead, are genuinely quantum fea- 
tures, beyond the quantum-to-classical correspondence, 
which may, however, maintain some tracts of universal- 
ity. Our analysis encompasses as examples the cases of 
quenches in the quantum Ising chain [7 and in a free 
bosonic system |15j. In addition, we show that, de- 
pending on the space dimensionality, the large deviation 
statistics of bosonic systems displays a so-called conden- 
sation transition (see, e.g., Ref. [16 ), analogous to the 
Bose-Einstein condensation. 

Statistics of the work - Consider a quantum system 
with N interacting degrees of freedom and Hamiltonian 
H{g). The extensive work Wm performed on the sys- 
tem during the quench go g is determined by the 
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FIG. 1: (a) Sketch of the excess free energy density fex{s) 
and (b) of the corresponding rate function I{w) for classes A 
(blue) and B (red) discussed in the main text. The gray area 
highlights the range of variables for which f^x does not have 
a thermodynamic interpretation. 



initial state, typically the ground state |^'q°) of H{go), 
and by the eigenvalues i?,^>o and eigenvectors j^'fj of 
the post-quench Hamiltonian H{g). In particular, Wn is 
a stochastic variable with probability density [iTt [T8] 

p{Wn) = E \{nm?mN - [ei - e^,°]), (i) 

n>0 

where indicates generically the extensive ground-state 
energy of H{g). As p vanishes identically for Wn below 
Ej^ - E^° , we refer Wn to this threshold so that Wn >0- 
The probability p(Wn) can be conveniently studied via 
its moment generating function 

G{s) ^ (e-^^"), (2) 

which, for N ^ oo, exists in the complex half-plane con- 
taining Res > (with possible zeros, see, e.g., Ref. 
For later purposes, we distinguish here a class A of sys- 
tems in which Wn for large but finite N cannot ex- 
ceed a certain extensive threshold Nwm from the class B 
whithin which Wn can assume arbitrarily large values. 
Generically, in the former class, G(s) is defined for all 
s e K with G(s) ~ e""^™^^ for s -oo, whereas in 
the latter G(s) is defined only for s > — s < with a 
generic singularity in its derivative at —s. The quantum 
Ising chain in a transverse field and the free bosonic field 
belong to classes A and B, respectively. 

The quantum to classical correspondence allows us to 
interpret the moment generating function G(s) for s > 
as the partition function of a classical system in a film 
geometry [9 . Indeed, Eq. ([ij implies 

G(s) = (*[J°|e-"[^(»)-^ol|4'g''), (3) 

where |e"''^(9) = Znxs is in fact such a par- 
tition function of the classical d + 1-dimensional system 
with transfer matrix e~^^^^ corresponding to the quan- 
tum Hamiltonian H(g), in a film geometry with trans- 
verse "surface" areaiV, "thickness" s and equal boundary 
conditions set by |^'o°). On the basis of Znxs one natu- 



rally defines the free energy Tnxs = — In-^TVxs per fcsT, 
where T is the temperature of the corresponding classi- 
cal system, which depends on the parameters of H{g). In 
terms of the classical system, the variable s in Eq. ([2| is 
the distance between the two confining surfaces which we 
assume to have a large transverse area N . Upon increas- 
ing s, the free energy density per unit area / decomposes 
in decreasing powers of s as |22] 

/ = N-^Fnxs ^ sfb + 2fs + corr., (4) 

where /{, = \vaiN.s^oo^Nxsl{Ns) is the bulk free en- 
ergy density and fs is the surface free energy density, 
i.e., the energy cost for introducing separately each sin- 
gle boundary into the otherwise bulk system. The cor- 
rections "corr." in Eq. ([4| vanish for s oo. In order 
to separate the effects of confinement from the bulk be- 
havior, one usually introduces the so-called excess free 
energy density per unit area f^x = f ~ sfb, which plays a 
fundamental role in what follows and becomes indepen- 
dent of N in the limit of large N considered hereafter. 
In terms of the quantum system, one finds from Eqs. ([3| 
and (|4|, that f, = E^/N, = -(In |vl/g) |)/iV 1 
and therefore 

G(s) =e-^^-(^l (5) 

For s < 0, fex is defined in terms of G(s) by this equations 
and it lacks its thermodynamic interpretation. 

Large deviations and universality - Equation ([s]) is 
crucial for understanding the emergence of universality 
in the large deviations statistics of the intensive work 
Wn = Wn/N. In fact, its distribution p{w) for oo 
can be determined by a saddle-point approximation of 
the inverse Laplace transform of G(s), which actually 
provides a heuristic derivation of the Gartner-Ellis theo- 
rem [2]. In particular, Eq. (|5| implies that p{w) has the 
form exp[—NI{w)], where the rate function I{w) is 
the Legendre-Fenchel transform of fex{s) (and viceversa, 
under certain assumptions [14j) 

I{w) = - inf {sw - fex{s)} , (6) 

in which the infimum is taken within the domain of def- 
inition of fex{s) and G{s). 

The generic features of p{w) can now be inferred from 
Eqs. (|5| and ([6|. First of all note that the excess free 
energy is such that fex{^) — and /ea.(0) = w. Most 
importantly fex{s) is a concave function of s |14| which 
approaches 2/^ for s — >■ +oo. Figure [T] provides a sketch 
of fex{s) and the corresponding /(s) for the two classes 
A and B introduced above. The last two properties im- 
ply the existence of a threshold in p{w): the infimum in 
Eq. ([6]) for It; < is —oo and consequently p{w < 0) = 0. 
The behavior of I{w) close to the threshold w > 0, in- 
stead, is determined by the one of fex{s) for s — !■ +oo 
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and in particular 7(0) = 2fs > 0, while the approach to 
it is determined by the corrections fex — 2fs in Eq. Q. 

The universality of these finite-size corrections close to 
critical points [21] carries over into the large deviation 
statistics of p(u') for w Ou. Indeed, if the post-quench 
Hamiltonian H(g) is close to a quantum critical point the 
finite-size corrections fex — '^fs to the free energy density 
of the (neax-critical) classical d + 1-dimensional system, 
which are responsible for the so-called critical Casimir 
effect |5T], take the universal scaling form s^'^0(s/^) for 
s ^ a, where ^ ^ a is the correlation length and a 
some microscopic length scale. The scaling function Q 
is universal in the sense of critical phenomena as 
it depends only on the universality class of the classi- 
cal critical point. In addition, due to the presence of 
the boundaries, Q depends on their surface universality 
class |23| or, equivalently, on which among the few ef- 
fective boundary (i.e., initial) states l^'o") flows 
to as the critical point is approached. Once the scaling 
function Q is known, the rate function is calculated via 
Eq. (|6|. In particular, if the post-quench Hamiltonian is 
critical, then = oo and 



liw < A) = 2/, 



1 



w \ 

a) 



d/(d+l) 



(7) 



with A = (i|9(0)|. While fs and the possible correc- 
tions depend on the specific parameters of the initial 
state, the leading dependence of I{w) on w is univer- 
sal and non-analytic. In the case of finite but large ^, 
the approach to the value 2/s is eventually controlled by 
Q{x ^ 1) = Cx°-e^^^ where the universal constants a, 
b, and C depend, along with Q, on the bulk and sur- 
face universality class of the transition, and they are 
known for a variety of universality classes [5T] (e.g., 
a = — 1/2,& = 2 for a quench of the quantum Ising 
chain within the same phase [9J). In this case, one finds 
I{w < ^-("i+i)) ~ 2fs ~ (C/6)wlnu;-i but with signif- 
icant logarithmic corrections. Notice that though this 
leading order is independent of a, a study of the rate 
function I{w) for larger values of w provides information 
on the complete scaling function Q, and hence on the 
corresponding boundary universality class. Accordingly, 
not only the edge singularities of the extensive work Wn 
studied in Ref. [9] are determined by universal features 
of the system and of the quench, but also the large de- 
viations of the intensive variable wn display universal 
properties in their rate function I{w) close to the thresh- 
old 

In order to illustrate the discussion above we focus on 
a free bosonic theory described by a Hamiltonian diago- 
nalizable in independent momentum modes 



H{m) 



where [(t>k,'!rk'] 



{2n) 



(8) 



first Brillouin zone \ki\ < tt. We assume a relati vis- 
tic dispersion relation ujk{m) = Vfc^ + rri^ and consider 
quenches of the mass from mo to m (H [15] . This Hamil- 
tonian captures the low-energy properties of a number 
of physical systems, including the ideal harmonic chain, 
interacting fermions and bosons in one dimension |26j . 
and it models the relative phase fluctuations of split 
one-dimensional condensates |3]. H{m) has a critical 
point at m = and the corresponding classical the- 
ory is that of a Gaussian field Lp m d + 1 spatial di- 
mensions and mass m. The quench is characterized 
by Afc = [wfe(mo) - LOk{m)]/[ujk{rno) + uik{m)] and from 
Eqs. ([T|, (U), and ([5| one finds [IS] 



fexis) = 



1 



d'^fc 

(27r)° 



In 



1 — ^2j^-2Ljfc(m)s 



(9) 



k,k' 



and the integral runs over the 



which is defined for s > — s = supj.(ln |Afe|)/a;fe(m) and, 
as anticipated, belongs to class B. This fex can be de- 
composed as in Eq. Q and upon approaching the critical 
point m = 0, i.e., for sufficiently large ^ — m^^ and s, the 
correction fexis) — 2/^ takes the (mo-independent) scal- 
ing form where Qo{^) is the scaling func- 
tion of the critical Casimir effect for the classical field 
if with boundaries belonging to the so-called ordinary 
surface universality class [23J, corresponding to Dirich- 
let boundary conditions for ip. Qq can be read, e.g., 
in Eq. (6.6) of Ref. [25]. Accordingly, upon approach- 
ing the critical point, the ground state Ivto"") of H{mo) 
flows towards the fixed-point state corresponding to 
this surface universality class. However, as s decreases, 
fex{s) — 2/s calculated from Eq. ([9| is no longer inde- 
pendent of mo and corrections to the critical Casimir 
term arise. These corrections are partly but effectively 
accounted for by changing s i— ^ s -|- 2i'oxt in the previous 
scaling form, where the so-called extrapolation length ^ext 
[23] takes here the value ttiq^. This can be interpreted 
as the fact that the fixed-point Dirichlet boundary con- 
dition on ip is effectively realized on confining surfaces 
which are at a distance ^oxt from the boundaries of the 
film of thickness s, resulting in an effective film of thick- 
ness s + 2^ext [H HI US] • For films sufficiently thick and 
critical, with s,^ ^ £oxt, this correction is unnecessary, 
while it becomes increasingly important as ^, s, or mo 
decrease. Figure [2] presents the rate function I{w) (solid 
line) in d = 1 for a quench from a non-critical to the 
critical point. In panel (a) the thin vertical dashed line 
indicates the mean work w. The dashed curve, instead, 
provides the quadratic approximation of I{w) around 
w = w, which describes the Gaussian distribution of the 
small fluctuations. While additional features of this I{w) 
are rationalized further below, we focus in panel (b) on 
the region of small w, where we expect universality to 
emerge. The dash-dotted line correspond to Eq. ([t]), with 
eo(0) = -r{d)Cid + l)/[{16Try/^r{d/2)] [25]. This uni- 
versal behavior sets in rather close to the threshold. How- 
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FIG. 2: Rate function I{'w) (solid line) of the work done on the 
lattice free bosonic theory in d = 1 and unit lattice spacing, 
for a quench from mo = 20 to the critical point. In panel (a) 
the dashed curve corresponds to the Gaussian distribution of 
small fluctuations around (w) — w. In panel (b), the dash- 
dotted curve is the prediction of Eq. ([7| while the dashed 
curve accounts also for a non- vanishing ^oxt. 



ever, the agreement between /(w) and Eq. (|7| extends to 
a wider range by accounting for the correction due to 
i'ext (dashed curve). The features displayed in Fig. [2] for 
m = carries over to the case m <^ niQ, which requires 
the knowledge of the full scaling function Qo{x)- For a 
fixed value of to, instead, the corrections to the scaling 
behavior due to ^cxt increase upon decreasing mg and 
eventually, after crossing the line to = mg of no quench, 
they lead to a change in the effective boundary state [5119] 
for mo —7- 0. 

Quantum regime and condensation - Let us now con- 
sider the case of large work w ^ w. Upon increasing 
w further away from the threshold, the value s*{w) of 
s for which the infimum in Eq. ^ is attained — and 
which satisfies /ea;(s*(w)) = w — decreases and so does 
the thickness of the corresponding film. The behavior of 
such a film is expected to become increasingly dominated 
by its microscopic details, with a generic lack of univer- 
sality even close to the critical point. Correspondingly / 
decreases because I'{w) — —s*{w). For w = w, s* = 
and I{w = w) vanishes, while it grows again for w > w, 
with s*{w) < (see Fig. [T|. 

The rate function for w > u" is thus determined by 
fex{s) for s < ("negative" film thickness), which lacks 
a thermodynamic interpretation because the quantum- 
to-classical correspondence does not hold in this case. 
The qualitative behavior of I{w > w) depends crucially 
on the class the system belongs to. In Fig. [l] we report a 
sketch of (a) fex{s) and (b) the associated rate function 
/ corresponding to classes A and B discussed above and 
characterized by (A) a bound (e.g., the quantum Ising 
chain) or (B) an unbound spectrum (e.g., free bosonic 
theory). In particular, in case A, I{w) diverges upon 
approaching wm, with I{w > wm) — +oo as required by 
the fact that p{w) vanishes above the intensive threshold 
Wm- In case B, instead, I{w — c») ~ sw and therefore 



p{'W ^ w) ^ e~^*™. This is seen in Fig. |2ja), though 
the asymptotic linear behavior for w ^ 1, with slope 
s (indicated by the dashed line) is actually approached 
only for rather large values of w. In general, s is system- 
specific and depends on the parameters of the quench. 

Even though the emergence of universality is appar- 
ently limited to w w, systems belonging to class B 
might display unexpected universal properties in the fully 
quantum regime w > w. In order to demonstrate this, we 
focus again on the free bosonic theory in Eq. (|8| and we 
show that for mg — ?> the statistics of the work displays 
a behavior analogous to the Bose-Einstein condensation 
of the ideal Bose gas in the grand canonical ensemble. 
This implies a transition in the large deviation statistics 
for w > w from exponential to algebraic. In fact, we 
note that the excess free energy fex{s) in Eq. ^ has the 
same form as half the scaled cumulant generating func- 
tion ip{s) of the fiuctuations of the spatial density pv of 
ideal Bose particles (of mass tob) within a large region of 
volume V . At equilibrium in an ensemble with chemical 
potential /i < (in units of temperature j3^^) one finds 

^(*) = / (CFln(i^A^), where A,. = 6-'^^^+^ with 

Ek = fi^k'^ /{2niB) and the integral is over E**. Accord- 
ingly the plot of ij{s) has the form B in Fig. [ija), with 
s = — /i. For the ideal Bose gas, the condensation occurs 
as /i — > /ic = 0: the asymptotic slope s — — /i of the rate 
function I{p > p) vanishes together with the function 
itself (see Fig. [ijb)). The mean value p = (p) = ^p'{0) 
above which this happens is the critical density for con- 
densation pc = l~'^(^{d/2) [57], which is finite only for 
d > dc = 2, where I = (27r/?fi^/TOB)^/^ is the thermal 
wavelength. / vanishes for p > pc because the probabil- 
ity p(p) acquires an algebraic dependence on p — due to 
the contributions of fiuctuations in single-particle states 
with small k — and indeed the momenta {{p~ p)") with 
n > d/2 diverge as /i — ^ pc', e.g., ((Ap)^) oc {—p)^" 
where a = 2 - d/2 for d < A. 

For the statistics of the work, TOg plays a role similar 
to p, although the occupation of the energy levels is de- 
termined by the non-thermal distribution generated by 
the quench and not by the Bose statistics. In fact, both 
TOg and p determine the fc-dependence of A^^q and A^^g, 
respectively, on which the onset of the condensation de- 
pends. In the case of the intensive work, TOq is the control 
parameter: for mg — >■ 0, (w) is finite for d > = 1 with 
a corresponding "critical value" Wc{m). The emergence 
of d'^ 7^ dc is due to the fact that the dependence of A^^g 
on k crosses over from quadratic for wg 7^ to linear for 
TOg ~ 0. Analogous crossover occurs in the condensation 
of an ideal Bose gas with relativistic dispersion [5S] . The 
rate function I{w) vanishes identically for w > Wc and 
p{w) acquires an algebraic dependence on w because of 
the slow asymptotic decay of the probability distribution 
of the work done on modes with small k, which are mildly 
confined in the initial state with TOg — >■ 0. As a result. 
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moments {{w — w)"') with n> d diverge in this hmit with, 
e.g., ((Aw)^) oc ln(m/mo) in d = 2. 

Conclusions - We discussed the qualitative features of 
the large deviation statistics of the work done during a 
quantum quench, highlighting the emergence of univer- 
sality and, for bosonic systems, of a non-thermal con- 
densation transition. Even though large fluctuations are 
exponentially rare as the system size increases, the value 
of the rate function /(O) can be reduced by a suitable 
choice of the quench parameters, making them observ- 
able by a post-selection of experimental data. 
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